Stochastic modeling of soil salinity

A minimalist stochastic model of primary soil salinity is proposed, in which the rate of soil salinization is determined by the balance between dry and wet salt deposition and the intermittent leaching events caused by rainfall events. The long term probability density functions of salt mass and concentration are found by reducing the coupled soil moisture and salt mass balance equation to a single stochastic differential equation driven by multiplicative Poisson noise. The novel analytical solutions provide insight on the interplay of the main soil, plant and climate parameters responsible for long-term soil salinization. In particular, they show the existence of two distinct regimes, one where the mean salt mass remains nearly constant (or decreases) with increasing rainfall frequency, and another where mean salt content increases markedly with increasing rainfall frequency. As a result, relatively small reductions of rainfall in drier climates may entail dramatic shifts in long-term soil salinization trends, with significant consequences e.g. for climate change impacts on rain-fed agricultur

An exactly solvable coarse-grained model for species diversity

We present novel analytical results about ecosystem species diversity that stem from a proposed coarse grained neutral model based on birth-death processes. The relevance of the problem lies in the urgency for understanding and synthesizing both theoretical results of ecological neutral theory and empirical evidence on species diversity preservation. Neutral model of biodiversity deals with ecosystems in the same trophic level where per-capita vital rates are assumed to be species-independent. Close-form analytical solutions for neutral theory are obtained within a coarse-grained model, where the only input is the species persistence time distribution. Our results pertain: the probability distribution function of the number of species in the ecosystem both in transient and stationary states; the n-points connected time correlation function; and the survival probability, definned as the distribution of time-spans to local extinction for a species randomly sampled from the community. Analytical predictions are also tested on empirical data from a estuarine fish ecosystem. We find that emerging properties of the ecosystem are very robust and do not depend on specific details of the model, with implications on biodiversity and conservation biology.Comment: 20 pages, 4 figures. To appear in Journal of Statistichal Mechanic

Pairing statistics and melting of random DNA oligomers: Finding your partner in superdiverse environments

Understanding of the pairing statistics in solutions populated by a large number of distinct solute species with mutual interactions is a challenging topic, relevant in modeling the complexity of real biological systems. Here we describe, both experimentally and theoretically, the formation of duplexes in a solution of random-sequence DNA (rsDNA) oligomers of length L = 8, 12, 20 nucleotides. rsDNA solutions are formed by 4L distinct molecular species, leading to a variety of pairing motifs that depend on sequence complementarity and range from strongly bound, fully paired defectless helices to weakly interacting mismatched duplexes. Experiments and theory coherently combine revealing a hybridization statistics characterized by a prevalence of partially defected duplexes, with a distribution of type and number of pairing errors that depends on temperature. We find that despite the enormous multitude of inter-strand interactions, defectless duplexes are formed, involving a fraction up to 15% of the rsDNA chains at the lowest temperatures. Experiments and theory are limited here to equilibrium conditions

Spatial effects on species persistence and implications for biodiversity

Natural ecosystems are characterized by striking diversity of form and functions and yet exhibit deep symmetries emerging across scales of space, time and organizational complexity. Species-area relationships and species-abundance distributions are examples of emerging patterns irrespective of the details of the underlying ecosystem functions. Here we present empirical and theoretical evidence for a new macroecological pattern related to the distributions of local species persistence times, defined as the timespans between local colonizations and extinctions in a given geographic region. Empirical distributions pertaining to two different taxa, breeding birds and herbaceous plants, analyzed in a new framework that accounts for the finiteness of the observational period, exhibit power-law scaling limited by a cut-off determined by the rate of emergence of new species. In spite of the differences between taxa and spatial scales of analysis, the scaling exponents are statistically indistinguishable from each other and significantly different from those predicted by existing models. We theoretically investigate how the scaling features depend on the structure of the spatial interaction network and show that the empirical scaling exponents are reproduced once a two-dimensional isotropic texture is used, regardless of the details of the ecological interactions. The framework developed here also allows to link the cut-off timescale with the spatial scale of analysis, and the persistence-time distribution to the species-area relationship. We conclude that the inherent coherence obtained between spatial and temporal macroecological patterns points at a seemingly general feature of the dynamical evolution of ecosystems.Comment: 5 pages, 5 figures. Supplementary materials avaliable on http://www.pnas.org/content/108/11/434

Towards a unified descriptive theory for spatial ecology: predicting biodiversity patterns across spatial scales

A key challenge for both ecological researchers and biodiversity managers is the measurement and prediction of species richness across spatial scales. Typically, biodiversity is assessed at fine scales (e.g. in quadrats or transects) for practical reasons, but often we are interested in coarser-scale (field, regional, global) diversity issues. Moreover, the pressures affecting biodiversity patterns are often scale specific, making multiscale assessment a crucial methodological priority. As species richness is not additive, it is difficult to translate from the scale of measurement to the scale(s) of interest. A number of methods have been proposed to tackle this problem, but most are too model specific or too rigid to allow general application. Here, we present a general framework (and a specific implementation of it) that allows such scale translations to be performed. Building on the intrinsic relationships among patterns of species richness, abundance and spatial turnover, we introduce a framework that links and predicts the profile of the species-area relationship and the species-abundance distributions across scales when a limited number of fine-scale scattered samples are available. Using the correlation in species' abundances between pairs of samples as a function of the distance between them, we are able to link the effects of aggregation, similarity decay, species richness and species abundances across scales. Our approach allows one to draw inferences about biodiversity scaling under very general assumptions pertaining to the nature of interactions, the geographical distributions of individuals and ecological processes. We demonstrate the accuracy of our predictions using data from two well-studied forest stands and also demonstrate the potential value of such methods by examining the effects of management on farmland insects across scales. The framework has important applications to biodiversity research and conservation practice

Emergence of structural and dynamical properties of ecological mutualistic networks

Mutualistic networks are formed when the interactions between two classes of species are mutually beneficial. They are important examples of cooperation shaped by evolution. Mutualism between animals and plants plays a key role in the organization of ecological communities. Such networks in ecology have generically evolved a nested architecture independent of species composition and latitude - specialists interact with proper subsets of the nodes with whom generalists interact. Despite sustained efforts to explain observed network structure on the basis of community-level stability or persistence, such correlative studies have reached minimal consensus. Here we demonstrate that nested interaction networks could emerge as a consequence of an optimization principle aimed at maximizing the species abundance in mutualistic communities. Using analytical and numerical approaches, we show that because of the mutualistic interactions, an increase in abundance of a given species results in a corresponding increase in the total number of individuals in the community, as also the nestedness of the interaction matrix. Indeed, the species abundances and the nestedness of the interaction matrix are correlated by an amount that depends on the strength of the mutualistic interactions. Nestedness and the observed spontaneous emergence of generalist and specialist species occur for several dynamical implementations of the variational principle under stationary conditions. Optimized networks, while remaining stable, tend to be less resilient than their counterparts with randomly assigned interactions. In particular, we analytically show that the abundance of the rarest species is directly linked to the resilience of the community. Our work provides a unifying framework for studying the emergent structural and dynamical properties of ecological mutualistic networks.Comment: 10 pages, 4 figure

Stochastic modeling of soil salinity

A minimalist stochastic model of primary soil salinity is proposed, in which the rate of soil salinization is determined by the balance between dry and wet salt deposition and the intermittent leaching events caused by rainfall events. The long term probability density functions of salt mass and concentration are found by reducing the coupled soil moisture and salt mass balance equation to a single stochastic differential equation driven by multiplicative Poisson noise. The novel analytical solutions provide insight on the interplay of the main soil, plant and climate parameters responsible for long‐term soil salinization. In particular, they show the existence of two distinct regimes, one where the mean salt mass remains nearly constant (or decreases) with increasing rainfall frequency, and another where mean salt content increases markedly with increasing rainfall frequency. As a result, relatively small reductions of rainfall in drier climates may entail dramatic shifts in longterm soil salinization trends, with significant consequences e.g. for climate change impacts on rain‐fed agricultur

Socio-Economic Instability and the Scaling of Energy Use with Population Size

The size of the human population is relevant to the development of a sustainable world, yet the forces setting growth or declines in the human population are poorly understood. Generally, population growth rates depend on whether new individuals compete for the same energy (leading to Malthusian or density-dependent growth) or help to generate new energy (leading to exponential and super-exponential growth). It has been hypothesized that exponential and super-exponential growth in humans has resulted from carrying capacity, which is in part determined by energy availability, keeping pace with or exceeding the rate of population growth. We evaluated the relationship between energy use and population size for countries with long records of both and the world as a whole to assess whether energy yields are consistent with the idea of an increasing carrying capacity. We find that on average energy use has indeed kept pace with population size over long time periods. We also show, however, that the energy-population scaling exponent plummets during, and its temporal variability increases preceding, periods of social, political, technological, and environmental change. We suggest that efforts to increase the reliability of future energy yields may be essential for stabilizing both population growth and the global socio-economic system

Statistical mechanics of ecological systems: Neutral theory and beyond

The simplest theories often have much merit and many limitations, and, in this vein, the value of neutral theory (NT) of biodiversity has been the subject of much debate over the past 15 years. NT was proposed at the turn of the century by Stephen Hubbell to explain several patterns observed in the organization of ecosystems. Among ecologists, it had a polarizing effect: There were a few ecologists who were enthusiastic, and there were a larger number who firmly opposed it. Physicists and mathematicians, instead, welcomed the theory with excitement. Indeed, NT spawned several theoretical studies that attempted to explain empirical data and predicted trends of quantities that had not yet been studied. While there are a few reviews of NT oriented toward ecologists, the goal here is to review the quantitative aspects of NT and its extensions for physicists who are interested in learning what NT is, what its successes are, and what important problems remain unresolved. Furthermore, this review could also be of interest to theoretical ecologists because many potentially interesting results are buried in the vast NT literature. It is proposed to make these more accessible by extracting them and presenting them in a logical fashion. The focus of this review is broader than NT: new, more recent approaches for studying ecological systems and how one might introduce realistic non-neutral models are also discussed

Association Patterns in Saproxylic Insect Networks in Three Iberian Mediterranean Woodlands and Their Resistance to Microhabitat Loss

The assessment of the relationship between species diversity, species interactions and environmental characteristics is indispensable for understanding network architecture and ecological distribution in complex networks. Saproxylic insect communities inhabiting tree hollow microhabitats within Mediterranean woodlands are highly dependent on woodland configuration and on microhabitat supply they harbor, so can be studied under the network analysis perspective. We assessed the differences in interacting patterns according to woodland site, and analysed the importance of functional species in modelling network architecture. We then evaluated their implications for saproxylic assemblages’ persistence, through simulations of three possible scenarios of loss of tree hollow microhabitat. Tree hollow-saproxylic insect networks per woodland site presented a significant nested pattern. Those woodlands with higher complexity of tree individuals and tree hollow microhabitats also housed higher species/interactions diversity and complexity of saproxylic networks, and exhibited a higher degree of nestedness, suggesting that a higher woodland complexity positively influences saproxylic diversity and interaction complexity, thus determining higher degree of nestedness. Moreover, the number of insects acting as key interconnectors (nodes falling into the core region, using core/periphery tests) was similar among woodland sites, but the species identity varied on each. Such differences in insect core composition among woodland sites suggest the functional role they depict at woodland scale. Tree hollows acting as core corresponded with large tree hollows near the ground and simultaneously housing various breeding microsites, whereas core insects were species mediating relevant ecological interactions within saproxylic communities, e.g. predation, competitive or facilitation interactions. Differences in network patterns and tree hollow characteristics among woodland sites clearly defined different sensitivity to microhabitat loss, and higher saproxylic diversity and woodland complexity showed positive relation with robustness. These results highlight that woodland complexity goes hand in hand with biotic and ecological complexity of saproxylic networks, and together exhibited positive effects on network robustness.The research Projects I+D CGL2011-23658 y CGL2012-31669 of the Spanish Minister of Science provided economic support